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Abstract — To monitor health information using wireless sensors 
on body is a promising new application. Human body acts as a 
transmission channel in wearable wireless devices, so electromag- 
netic propagation modeling is well thought-out for transmission 
channel in Wireless Body Area Sensor Network (WBASN). In 
this paper we have presented the wave propagation in WBASN 
which is modeled as point source (Antenna), close to the arm 
of the human body. Four possible cases are presented, where 
transmitter and receiver are inside or outside of the body. Dyadic 
Green's function is specifically used to propose a channel model 
for arm motion of human body model. This function is expanded 
in terms of vector wave function and scattering superposition 
principle. This paper describes the analytical derivation of the 
spherical electric field distribution model and the simulation of 
those derivations. 

Index Terms — Wireless Body Area Networks, Dyadic Green's 
Function 



I. Introduction 

Hospitals throughout the world are facing a unique problem, 
as the aged population is increased, health-care population is 
decreased. Telecommunication community is not doing much 
work in the field of medicine however, there is a need of 
remote patient monitoring technology. To fulfill this task, 
it is required to build communication network between an 
external interface and portable sensor devices worn on and 
implemented within the body of the user which can be done 
by BASNs. 

BASNs is not only useful for remote patient monitoring, 
but can also establishes within the hospitals; like in operation 
theaters and intensive care units. It would enhance patient 
comfort as well as provide ease to doctors and nurses to 
perform their work efficiently. BAN is used for connecting 
body to wireless devices and finds applications in various areas 
such as entertainment, defense forces and sports. 

The basic step in building any wireless device is to study 
the transmission channel and to model it accurately. Channel 
modeling is a technique that has been initiated by a group of 
researchers throughout the world [1]. They have studied path 
loss and performed measurement campaigns for wireless node 
on the body [2-8]. Some researchers have taken into account, 
the implanted devices which are the area of BAN called 
as intra-body communication [9]. For the short range low 



data rate communication in BAN, measurement groups have 
considered Ultra-Wide Band (UWB) as the appropriate air 
interface. The models developed by measurement campaigns 
are only path loss models and do not provide any description 
of propagation channel. 

It is important to study the propagation mechanism of radio 
waves on and inside the body in order to develop an accurate 
BAN channel model. This study will show the underlying 
propagation characteristics. It would help in the development 
of BAN transceivers which are much suited to the body 
environment. 

For a given position of the transmitter on or inside the body 
it is required to find out the electromagnetic field on or inside 
the body for a BAN channel model. This is quite a critical 
problem that requires a large amount of computational power. 
Therefore, it is necessary to derive an analytical expression 
which will perform this objective. In short this determines 
which propagation mechanism takes place, that is reflection, 
diffraction and transmission [10]. An appropriate method of 
doing this task is by using Dyadic Green's function. The 
solution of canonical problems, such as cylinder, multi layer 
and sphere have been solved in Electro Magnetic (EM) theory, 
using Dyadic Green's Functions [11-13]. 

II. Motivation 

Recently, WBASNs shows potential due to increasing ap- 
plication in medical health care. In WBASNs, each sensor 
in the body sends it's data to antenna,both sensors and 
antenna are worn directly on the body. Examples include 
sensors which can measure Brain activity, blood pressure, 
body movement and automatic emergency calls. We require 
simple and generic body area propagation models to develop 
efficient and low power radio systems near the human body. To 
achieve better performance and reliability, wave propagation 
needs to be modeled correctly. Few studies have focused on 
analytic model of propagation around a cylinder (as human 
body resembles a cylinder) using different functions. These 
functions involve Mathieu function, Dyadic Green's function, 
Maxwell's equations, Finite Difference Time Domain (FDTD) 
and Uniform Theory of Diffraction (UTD). Some of these 
approaches have already proven effective for evaluating body 
area communication system proposals. 
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Finite Difference Time Domain had successfully measured 
the communication scenarios. Complete Ultra-Wide band 
models have been developed using measurements and simula- 
tions, however they do not consider the physical propagation 
mechanism. So, the researchers have to rely on ad-hoc model- 
ing approaches which can result in less accurate propagation 
trends and inappropriate modeling choices [14, 15]. 

Uniform Theory of Diffraction depends on a ray tracing 
mechanism allowing propagation channel to be explained in 
terms of ray diffraction around the body . It typically based 
on high-frequency approximations which is not valid for low 
frequencies, also not useful when antenna is very close to the 
body [16]. 

A generic approach is proposed to understand the body area 
propagation by considering the body as a lossy cylinder and 
antenna as a point source by using Maxwell's equation. A 
solution for a line source near lossy cylinder is derived using 
addition theorem of Hankel functions then the line source 
is converted into the point source by taking inverse Fourier 
transform. The model accurately predicts the path loss model 
and can be extended to all frequencies and polarities but this 
is limited in scope and not always physically motivated [17]. 

Mathieu functions are also used for body area propagation 
model. The human body is treated as a lossy dielectric elliptic 
cylinder with infinite length and a small antenna is treated as 
three-dimensional (3-D) polarized point source. First the three- 
dimensional problem of cylinder is resolved into 2-D problem 
by using Fourier transform and then this can be expanded in 
terms of Eigen functions in cylindrical coordinates. By using 
Mathieu function exact expression of electric field distribution 
near the human body is deduced [18]. 

The propagation characteristics of cylindrical shaped human 
body have been derived using Dyadic Green's functions. 
The model includes the cases of transmitter and receiver 
presents either inside or outside of the body and also provides 
simulation plots of Electric field with different values of 
angle (6). All the above proposals describe the propagation 
characteristics of cylindrically shaped human model [19]. 

We have developed a simple but generic approach to 
body area propagation derived from Dyadic Green's Function 
(DGF). This approach is for arm motion of human body. 
When the human arm is moved in r, 6, (f> direction, propagation 
characteristics of spherical shaped have been derived using 
DGF. First, we use spherical vector Eigen functions for finding 
the scattering superposition. Four cases are considered for 
either transmitter or receiver is located inside or outside the 
body. Finally, simulated results of electric field distribution 
with different values of angle have shown. 

III. Mathematical Modeling for Arm Motion 
using Dyadic Green's Function 

In this paper, spherical symmetry is used to represent in 
and around the arm of the human body. A point on body is 
a sensor, denoted by x which represents (r,Q,(j>) coordinates 
in the spherical coordinate system and xo is the location of 
transmitting antenna. (r,0,</>) are unit vectors along radial, 
angle of elevation from z-axis and azimuthal angle from x- 
axis as shown in figure 1. 




x 



Fig. 1. Human body model showing arm motion in 3D. 

A. Electric Field Propagation Characteristics 

Let E{x) be electric field at point x due to current source 
J(xq). The general formula for Electric field can be written 

as: 

E(x) = iuifip J J J G(x,x )J(x,xo)dv (1) 

V is volume of source, J(Xq) is the current source, 
G(x,xo) is the Dyadic Green's function 'u/ is the radian 
frequency of transmission and ' ' p! p is magmatic permeability of 
the medium. A Dyadic Green's function is a type of function 
used to solve inhomogeneous differential equations subject to 
specific initial conditions or boundary condition. 

B. Spherical Wave Vector Eigen Function 

As we are considering arm motion of human body, so 
spherical symmetry is used by taking shoulder as center. For 
this, spherical eigen functions are used to write the Dyadic 
Green's function. 

Dyadic Green's function is basically depends on the spher- 
ical vector eigen functions [14]. These eigen functions are 
L n hk(x)> Mnhk(x) and N nhk (x), where k is the wave number 
of medium, n is an integer, h is a real number and a; is a point 
in space. These all are the solutions to the Helmholtz equation 
having three components in r, 9 and <fi. These vector eigen 
functions are given by [19]: 

Lnhkix) = V[* nh fc(x)] (2) 



M nhk { X ) = V x lnhkix)} (3) 



N nhk {x) = y;V x V[tf nWfc (x)] (4) 
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In above eigen functions, Laplacian operator in the spherical 
coordinate system is V. It's mathematical expression is given 
as: 



d_ 

dr 



d_ 

7d0 



d 



r sin 



(5) 



x represents the point in space having components r, Q 
and <f>. Solution of Helmoltz equation is ^ n hk{x) which is the 
scalar eigen function [19]. 



[*nfcfc(x)] = Z n ( V r)P% (cos 6)Zh<t> 



(6) 



Z n is a general spherical function of order n. For sphere 
we use Hankie function of first and second order which are 
defined as: 



C. Scattering Superposition 

In scattering problems, it is desirable to determine an 
unknown scattered field that is due to a known incident field. 
Using the principle of scattering superposition we can write 
Dyadic Green's equation as superposition of direct wave and 
scattering wave. In Figure 2, concept of scattering superposi- 
tion is shown in which there is a sensor located inside the arm 
of body considered as sphere. The sensor transmits the wave 
to antenna which is divided in two parts as Direct wave and 
Scattered wave. The Direct wave is considered as wave directly 
transmits from sensor to transmitter and scattered wave is 
composed of reflection and transmission waves. Therefore, 
general equation of scattering superposition is illustrated as: 



G(x, x ) = G d (x, x ) + G s (x, x ) 



(14) 



rj is the propagation constant in direction of <fi, whereas 
k 2 = rj 2 + h 2 . The laplace operator is applied and find the 
eigen values L nhk , M nhk and N nhk by using eigen function. 
The vector eigen function in (2), (3) and (4) becomes: 



(8) 



M nhk {x) 



T hZ^ p h {coa0) ^_ Zn{vr) 



sin 6 



89 



P n h (cos0)» 



(9) 



N nhk ( x ) = ^^P^cosemH + 1 



(10) 



These three vector eigen function are perpendicular among 
themselves as well as with respect to each other [11]. In the 
form of matrices, vector Eigen functions can be written in this 
form, 



L n hk{x) 



M nhk { X ) 



h ^Mph( C0S ey, 



sin 




I 



o 



V -Z n {w) 



dP* (cos 8)lZh4> 

de 



(11) 



(12) 



Dyadic Green's equation is divided in to two parts as direct 
wave [Gd(x,x )] and scattered wave [G s (x, Xo)]. The direct 
wave corresponds to direct from source to measuring point 
and scattered is the reflection and transmission waves due to 
presence of dielectric interface. 



D. Superposition of Direct Wave 

The direct component of DGF is given as [11]: 



G d (x,x ) ■■ 



rr 



+ 1 



dh 



\M^ k {X)®M* hk {X ) + NV k { 
{ M nhk (X) M^*l (X ) + N nhk ( 



-oo 



n— — oo 



n 2 



-xx 



■.(Xo)+N^ k (X)®N: hk (X a ) 

■AX)®N^ k \x ) 



(15) 



In the above equation of DGF, r > is for first case and 
r < ro is second case. The * denotes the conjugation and (g) 
is for the Dyadic product. Here we introduces superscript (1) 
for outgoings wave and other for standing waves. If the vector 
eigen function has the superscript (1) then, is chosen for 
Z n and J n should be used otherwise. 



E. Superposition of Scattered Wave 

Here we discuss four different scenarios for the scatter- 
ing components of DGF along with boundary conditions 
G s (x,x ). (i) Both receiver and transmitter are inside the 
body, (ii) The receiver is located outside and transmitter is 
located inside the body, (iii) The receiver is located inside 
and transmitter is outside the body, (iv) Both transmitter and 
receiver are located outside the body. Receiver and transmitter 
are in the order: 1 denotes the medium inside human body 
and 2 is for free space medium. 



Nnhkix) 



d n {r]r) P% {cos 0)1?* h<f> 
krdr d0 



hcj) 



(13) 



F. Transmitter and Receiver Located Inside Body 

In this case, Receiver and Transmitter both located inside 
the body so we can write Dyadic Green's equation as, 
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Fig. 2. Scattering Superposition 



J -°° -oo ' 
\ht i\ t Inn \ N nhkl (X ) T 

x[M nhk iN nhk i]Ri2 x i 



(16) 



i?i2 contains reflection coefficients. R\2 is calculated in 
literature using boundary conditions, its matrix is given by 
[16]: 



In the above equation T\2 is a transmission coefficient 
Matrix and given as: 



T12 



2uj 



[J n {7]-i_d)H n (7] 2 d) - H n {7] 2 d)J n {r]id)\ 1 



£1 



£ 



(20) 



H. Both Transmitter and Receiver Located Outside Body 



(17) 



R12 = [J n (vid)H n (r]2d) - H n (r)2d)J n (r)id)]~l 

x[H n {q 2 d)H n {i ll d) - if„(r?id) J n (r} 2 d)]~l 

In the above equation of reflection coefficient 'd' represents 
radius of spherical body model, rjf = h\ — h? = k\ — 
h 2 ,kf = w 2 /zi£i,fc| = uj 2 fi2^2- The 2x2 matrices for j n {r]d) 
and H n (rjd) are expressed as: 



B n (rjpd) 



1 ( jue p 7] p dB n (ripd) 
r&d V -nhB n (r] p 



-nhB n (rj p 
-juj[i p ripdB n (r)pd) 



G s (x,x ) 



/OO 00 1 



n— — 00 



XlMnhfeJVnhfcJ^l 

{iV„/ l A:(Xo) T M n / l A : (Xo) T 



(21) 



B n is either Hn° rJ "\ B(.) is the derivative of B w.r.t the 
whole argument, and p=l,2 

G. Transmitter Located Inside and Receiver Located Outside 
Body 

In this case DGF can be written as : 



Similarly as R12, R21 is the reflection coefficient matrix 
and it is given as: 
(18) 

i?21 = [J n { Vl d)H n { m d) - H n ( V 2d)J n (riid)}-l 



x[Jn(.md)Jn{riid) - J n (r]id)Jn(r]2d)] 



(22) 



/. Transmitter Located Outside and Receiver Inside Body 
In this case, we can write DGF as: 



/OO 00 -f 

^ 71 — — OO 

x[N nhk .M nhk ]T 12 (^i hkl t°Z 

\ M nhkl{ x o) 



/OO OO ^ 



(19) 



x[M nhk N nhk }T 21 

^ M nhkl( X of 



(23) 
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T12 is the transmission coefficient matrix, given as: 



T21 



2cj 



[Jn(m d )H n (f] 2 d) - H n (r] 2 d)J n (r]ld)] 1 



£2 






-M2 



(24) 



IV. Transmitter and Receiver Located Outside of 
the Body 

In this section we presents the equation which is required 
for simulation. With the help of simulation it will be easy to 
study the propagation characteristics of arm motion making 
spherical pattern. 



G s 



dh 



00 ^ 

V" -^G nh (x,x )dh 

n— — 00 

G n h(x 7 x ) is stated as: 



(25) 



L and Q are the truncation limits and AH are the step 
size of integration. N and AH are so small that could be 
ignored and has no effect on calculations. We only presents 
electric propagation of multi-path reflection and transmission 
waves of scattering DGF.This is more significant to represent 
the attribute of arm motion as compared to the direct DGF. 
Figure 2,3 and 4 show the scattering DGF (simulation) of 
electric field with the change in 9. 

phi vs magnitude vs distance 



10 



n 6 



4 




G nh (x,x ) = (NnhkiX) 1 M nhk {Xf) x R21 

(N nhk (XofM nhk (X ) T ) { ' ib) 

V. Simulations 

As we have defined earlier, arm motion at different angles 
are presenting spherical pattern. Therefore, we simulate the 
radio propagation environment having radius d = 15cm, meg- 
natic permeability for human body (assume that permeability 
of human body is approximately equal to air) fi 2 = 1-256 x 
10~6, similarly electric permittivity e 2 = 2.563 x 10~10. 
The dielectric constant is mean value of all tissues of human 
body. We take the surrounding homogeneous medium to be 
air with megnatic permeability [i\ = 1.256 x 10~6 and 
electric permittivity e\ = 8.8542 x 10~12. Frequency up 
to GHz is used for BAN communication, which is for ISM 
band. The Transmission frequency for simulation is 1GHz. 
We assumed that the transmitter is acting as point source 
at x a = (16cm, §,0). The radial distance of receiver is 
r = 18cm from the central spherical axis of shoulder. For the 
simulation, we assumed that receiver move along the azimuthal 
angle for varying values of O and different heights from the 
center of shoulder. 

For simulation, we consider equation (25) in which 
G n h(x,x ) is used in matrix form of eigen functions. This 
equation has an integration which is not possible so we 
approximate it to summation. Thus, we approximate equation 
(25) in to this form: 

L 

Gs(x,x a ) = — 
V* -^G nh {x,x )dh 

n=-Q 




P hl ( radlans ) 4Ubb magmtude(dB) 

Fig. 3. Magnitude of scattered field component versus angle 0,with 
different values of d and the angle is 6 = ^ 

Using equation (27), we have three components in r,9 and 
(f> direction. Every Component of electric field is plotted as a 
function of azimuthal angle (f>. The values of is (0 to 2n), 
whereas at z coordinate different values of receiver has been 
plotted. The electric field is plotted, which is vector addition 
of three components. These all parameters are shown in the 
simulation graph. 

By taking the value of — |, figure 2 shows that magnitude 
of electric field (E^) is decreasing as the distance of receiving 
antenna is increasing from the sensor (transmitting antenna). 
The plot shows electric field component at different values of 
</>, varying from to 2-k. In this case, E$ is decreasing from 
(4080 to 4065)dB by replacing the receiving antenna from 
cm to 10 cm. 

In Figure 3, when we take value of 9 = f, magnitude of 
electric field (E^) again decreases as the antenna moves away 
from sensor. For the values of <p from to 27r, E^ has different 
values from (4060 from 4068)dB. By changing position of 
receiving antenna from cm to 10 cm. 

The values of distance and cf) are same, as described in the 
above graphs by only replacing the parameter 9 = n. Similarly 
in figure 4 values of E^ change from (4082 to 4074)dB 
by moving the position of receiver away from transmitting 
antenna, which in return decreases the electric field intensity. 
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magnitude vs phi vs distance 



10 v 




distance(cm) ™ magnnltude (dB) 



Fig. 4. Magnitude of scattered field component versus angle <^>,with 
different values of d and the angle is 9 = ^ 



magnitude vs phi vs distance 



10 v 




distance(cm) «4 magnnltu de(dB) 

Fig. 5. Magnitude of scattered field component Em versus angle <^>,with 
different values of d and the angle is 8 = ir 

VI. Conclusion 

We have proposed a generic approach to derive an ana- 
lytical channel modeling and propagation characteristics of 
arm motion as spherical model. To predict the electric field 
around body, we have formulated a two step procedure based 
on Dyadic Green's function. First, we derive Eigen functions 
of spherical model then calculated the scattering superposition 



to come across reflection and transmission waves of antenna. 
The model includes four cases where transmitter or receiver is 
located inside or outside of the body. This model is presented 
to understand complex problem of wave propagation in and 
around arm of human body. Simulation shows that Electric 
field decreases when receiver moves away from the shoulder 
with change of angle 0. 
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